Renormalization group, dimer-dimer scattering, and three-body forces 
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We study the ratio between the fermion-fermion scattering length and the dimer-dimer scattering 
length for systems of nonrelativistic fermions, using the same functional renormalisation technique 
as previously applied to fermionic matter. We find a strong dependence on the cutoff function used 
in the renormalisation flow for a two-body truncation of the action. Adding a simple three-body 
term substantially reduces this dependence. 
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Ultra-cold Fermi gases provide a fertile ground for re- 
search in atomic physics. An important feature of such 
systems is superfluidity, which is the result of attractive 
fermion-fermion interactions leading to pairing. Recent 
advances using Fcshbach resonances allow a tuning of the 
fermion-fermion S-wave scattering length ap. For neg- 
ative scattering length we get the weak-coupling BCS 
state. For positive values of ap bound states of two 
fermions- "dimers" -form and these can lead to a Bosc- 
Einstein condensate (BEC) [l|. The size of dimers is 
determined by the fermion-fermion scattering length and 
their binding energy is of order l/a 2 F . 

For a sufficiently dilute and cold gas of dimers the 
main dynamical quantity characterising their interaction 
is the dimer-dimer scattering length as- The exact rela- 
tion between dimer-dimer and fermion-fermion scattering 
lengths a,B = 0-6ap was established in Ref. [j| by solving 
the Schrodinger equation for two composite bosons in- 
teracting with an attractive zero-range potential. This 
method is difficult to extend to the many-body case. 
Therefore, it is useful to study the ratio flg/ajf in an ap- 
proach which can be used both for few and many-body 
problems. 

In this paper we calculate as in the framework of the 
Exact Renormalisation Group (ERG) approach. (For re- 
views, see Refs. [1, H|). This technique has been previ- 
ously used to study a variety of physical systems, from 
systems of nonrelativistic fermio ns fB -H, [ToL [Tl} to quark 
models [l2| and gauge theories [i~3l]7 It is based on the 
scale-dependent average effective action , where k is an 
auxiliary running scale. The action at scale k contains 
the effects of field fluctuations with momenta q larger 
than k only. In the limit k — > all fluctuations are in- 
cluded and the full effective action is recovered. In prac- 
tice one introduces a set of fc-dependent cutoff functions 
R(q), which suppress the effect of modes with q < k in 
the path integral for the action by giving them a large 
k dependent mass. The functions R(q) should vanish in 
the limit k — > and behave like k a with a > when 
q -> 0. 

With this prescription the average effective action at 
large k is just the classical action of the theory — in our 
case nonrelativistic fermions with a local interaction. An 



exact solution of the functional RG equation should be 
independent of the choice of cutoff for k — s- 0. However, 
in practice, truncations of the action inevitably lead to 
some cutoff dependence of the results. We can use this 
dependence as a measure of the quality of the truncation. 
With this tool, we shall see that the standard parametri- 
sation of the effective action containing only two-body 
terms is insufficient, and that we get better results by 
including the simplest three-body term. 
The flow of the effective action satisfies 131 
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(d k R) (r( 2 > - R)- 
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where T*- 2 ' is the second functional derivative with re- 
spect to the fields, and the cutoff functions in the mass- 
like term R(k) drive the RG evolution. The opera- 
tion STr denotes the supertrace [13] taken over energy- 
momentum variables and internal indices and is defined 
by 

The evolution equation for the average effective action 
has a one-loop structure, but contains a fully dressed, 
scale-dependent propagator (r^ 2 ) — R)~ x . Thus, despite 
its apparently simple form, Eq. ([I} is actually a func- 
tional differential equation. In the absence of general 
methods to solve such equations numerically we must 
resort to approximations. One common approach is to 
parametrise the effective action with a finite set of terms, 
turning the evolution into a system of coupled ordinary 
differential equations for their coefficients. These equa- 
tions can then be solved numerically. Here we study pos- 
sible truncations for fermionic few-body systems, and our 
choice of ansatz for the action is motivated by both ERG 
studies of many-body systems 0, @ and effective field 



theories (EFTs) for few-fermion systems [15j. The tech- 
nique we use is similar to the one used in Ref. 0] to 
analyse the scattering in the system of two nonrelativis- 
tic fermions. 

A rather different approach to the low-energy fermion- 
dimer scattering was considered in Ref. [8|. There 
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the Skomiakov-Ter-Matirosian equation 0, [Tsj] was de- 
rived from the RG flow with energy- and momentum- 
dependent three-body couplings. Whilst the results ob- 
tained in that work demonstrate the formal equivalence 
between an RG approach and three-body quantum me- 
chanics, it seems to be very difficult to extend the treat- 
ment to the more complicated cases of four-body or 
many-body systems. In the present paper we focus on 
the effect of the three-body forces in the system of two 
dimers, keeping in mind possible extensions of the for- 
malism to many-body systems, such as that attempted 
in Ref. . Therefore we only include the simplest three- 
body term. 

The formation of the trion - the correlated state of 
three fermions - was also studied in Rcfs. dd, d, but 
the effect of the three-body interactions in the context of 
the four-body systems was not considered. 

We first examine the case when fermions interact 
only pairwise. This has previously been considered in 
Refs. @, Q, where many the technical details can be 
found, and so we give only a brief account of the formal- 
ism, concentrating on the dimer-dimer scattering length 
a b ■ We use this to extend the results of Ref. [(| and ex- 
amine the cutoff dependence of as- We then turn to our 
main task, the inclusion of three-body terms in T. 

As the cutoff scale tends to infinity, we demand the 
action to be a purely fermionic one containing a contact 
two-body interaction without derivatives. This kind of 
interaction has been extensively used in the EFT-bascd 
studies of nuclear forces [l5[. It is convenient to rc- 
cxprcss the theory in terms of an auxiliary composite 
boson field by making a Hubbard-Stratonovich trans- 
formation. This replaces the two-body interaction by 
a Yukawa-type coupling between the fermions and the 
auxiliary boson. A kinetic term for the boson is then 
generated by the RG evolution. The minimal effective 
action used in previous work is 0- 
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Here H(x,x',k) is the scale-dependent boson self-energy 
and 1*2 parametrises the boson-boson interaction which 
can be generated by the evolution. The latter is equiva- 
lent to a four-body interaction in terms of the underlying 
fermions. To this action we add a local three-body inter- 
action, similar to that used in another context in Ref. fioj 
Expressed in terms of the boson field this has the form 



-A / d A x^(x)^(x)cf>(x)ip(x) 



(4) 

We concentrate first on the two-body part of Eq. ([3]) . 
The evolution of the boson self-energy is given by 



d k n( x ,x',k) 



S(j){x')S^ (x) 



(5) 



although from now on we shall express all evolution in 
momentum space. Note that only fcrmion loops con- 
tribute to the evolution of the boson self-energy in vac- 
uum. These depend on the fermionic cutoff function Rp, 
for which we take the form I17l 



R F {q,k) 



q 



2M 



6{k-q). 



(0) 



We impose the boundary condition that the scattering 
amplitude in the physical limit k — > reproduces the 
fcrmion-fcrmion scattering length, 



(7) 



Here Pq (P) denote the total energy (momentum) flow- 
ing through the system and p = y / 2MPo — P 2 /2 is the 
relative momentum of the two fermions. Integrating the 
resulting ERG equation gives Q 
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Using a gradient expansion of the action, we can define 
boson wave-function and mass renormalisation factors by 



and 
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where So = — l/(^ a F) denotes the bound-state energy 
of a pair of fermions. Note that these renormalisation 
factors are only identical in vacuum for a limited sub- 
set of cutoff functions (which, like (JSJ) must preserve 
Galilean invariance to lowest order), otherwise the iden- 
tity Z^(k) = Z m (k) holds only in the physical limit k — >• 
and the evolution of these renormalisation factors should 
be calculated separately. 

The evolution of the boson-boson scattering amplitude 
follows from 
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-d k u 2 (£D,k) 



S^(£ D ,0)6^(£ D ,0) 



<9aT|0 = o. 
(11) 



This equation can be separated into fcrmionic and 
bosonic contributions containing dkRF an d d k RB, re- 
spectively. We first look at the mean-field result, where 
bosonic contributions are neglected. The evolution of u 2 
is then given by 



dku 2 



3^ 
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d'-t 



dkRF 
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where Epn(q, k) 
gives 



(2tt) 3 [(E FR (q,k)-£ D /2)Y 
2jj q 2 + Rp{q,k). Integrating this 



u 2 (o) = — mV4. 
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where we have again used the sharp cutoff function of 
Eq. ([6]). The scattering amplitude at threshold is 



Tbb = 



8tt 2w 2 (0) 
a B = 
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giving the well-known mean-field result 0,3 = 2clf [18j 
which is far from the exact value as = 0.6cif @- This im- 
plies that beyond-mean-field effects such as dimer-dimcr 
rescattcring arc important be considered. 

To include such effects we must take into account the 
boson loops. After some algebra, we find 



dku 2 \ 
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d 3 t 



d k R B 



2Z 3 (k)J (2tt) 3 [E BR (q,k)-£ D f 



where 



E BR (q,k) = —q 2 



ui(k) R B {q,k) 
Z<f,(k) Z${k) 



and 



iti(fc) = -U(£ D ,0,k). 



(15) 
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We choose the bosonic cutoff function to be as close as 
possible to the fermionic one, 



R B (q,k) = Z (i 



{c B kf - q 2 
AM 



(c B k - q), 



(18) 



apart from the addition of a parameter c B , which sets the 
relative scale of the fcrmionic and bosonic regulators, and 
a factor of Z^. The latter has the important advantage 
of leading to a consistent scaling behaviour, so that all 
contributions to a single evolution equation decay with 
the same power of k for large k. Moreover it also gives 
cif-scaling, where all terms in a single equation have the 
same dependence on aj?. 

The mean-field result is recovered for c B = while 
the opposite limit c B — > 00 leads to a B — > 0. Using 
c B = 1 gives a ratio of a B /a,F = 1.13. Taking c B = v2 
as in Ref. @ results in a B /ap = 0.75. In general, the 
results show a rather strong dependence on the relative 
scale parameter c B , as can be seen below in Fig. [5J Such 
dependence of a physical result on the choice of cutoff 




FIG. 1: The skeletons of the diagrams that contribute to the 
evolution of the three-body term. A dashed line denotes a 
boson, a solid line a fermion. Each diagram can have one 
insertion of dkRF.B as appropriate on a single internal line. 
We have the ladder diagram on the left, two triangle diagrams 
and finally the box diagram. With insertions there are 24 
diagrams in total. 



must be an artifact of our truncation of the action, since 
it should vanish for an exact solution of the functional 
RG equation. It indicates that the minimal ansatz used 
for the effective action needs to be extended. 

We now consider the effect of adding the local three- 
body force of Eq. (|4]). The evolution equations for iti and 
Z<i, remain unchanged but the one for u 2 now becomes 



dkU 2 = t 
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The evolution equation for A is defined by an expansion 
about the energy of the bound state pole for bosons, and 
half that energy for fermions, 



d k X = -- 



<5 4 STr [d k R(rW _ R)~ 



2 64>XS D ,0)S4>(S D , 0)Si/>\e D /2, 0)6ip(£ D /2, 0) ' 

(20) 

There are three distinct contributions to the running 
of A, coming from ladder, triangle and box diagrams, as 
shown in Fig. [T] We denote the corresponding driving 
terms as Di, D t and Db, splitting the last two into their 
fermionic and bosonic contributions. After evaluation of 
traces and contour integrals, we get 



Di 



D t 



D 



A 2 



9 2 X 



9 2 X 



d 3 q dk(R F Z^+dkR B 



(2tt) 3 {Ef^dZ^ + Eb^dY' 
d 3 q dk,R F {E BRi D + 2Z < i > Efr,d) 
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d 3 q 
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dkRs 
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Ebr,d) 2 



(21) 
(22) 
(23) 
(24) 
(25) 
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FIG. 2: Ratio of boson-boson to fermion-fermion scattering 
lengths as a function of the relative scale parameter cb. The 
blue (solid) curve shows results for the minimal action r m i n ; 
the red (dashed) curve shows the effect of adding the local 
three-body term. 

where E F r,d = E F r - £ D /2 and E B r.d = E B r - Ed- 
Note that the evolution of the three-fermion interaction 
A does not depend on the four- fermion interaction u 2l 
but u 2 does depend on A. For the initial condition on A, 
we demand that there is no fundamental three-fermion 
interaction, and hence A — > as k — > 00. In practical 
calculations it is simpler to take A to be zero at some 
large starting scale. Since A behaves like 1/fc 2 for large 
k, this can be used reliably. 

Now we turn to the results. We note that these are 
numerically independent of the starting scale provided 
it is chosen to be at least ka F — 100. For k ^> 1/ap 
the system is in the "scaling regime" , and a fixed point 
corresponding to the unitary limit governs the evolution 
until k becomes comparable with 1/ap- (This can best 
be seen from the evolution of the dimensionless coupling 
k 2 \.) 

Wc find that the ratio as jap decreases when the three- 
body term is included. For example, choosing eg = 1 
leads to as jap = 0.74, which should be compared to 
clb/qf = 1-13 without the three-body term. Similarly 
the choice cb = V% gives asjap — 0.69. 

The full behaviour of as /ap as a function of cb is pre- 
sented in Fig. [2j This shows that, as well as reducing 
the overall size of the ratio (Ib/clf, the inclusion of the 
three-body force significantly weakens its dependence on 
the relative scale cb- We expect the qualitative features 
of this picture to remain correct for any bosonic regula- 
tor although the quantitative details will depend on the 
particular functional form used. 

Note that for large cb, the dominant contributions to 
clb/clf comes from the boson-loop terms in the equa- 
tion for u 2 . Since these do not depend on the three- 
body coupling A, the two curves approach each other. 
Moreover, this limit corresponds to integrating out the 
fermions first, which generates a non-zero value for u 2 at 



the start of the bosonic integration. In the limit cb — > 00, 
this coupling is driven to the trivial fixed point, u 2 = 0, 
since we have no terms to cancel the linearly divergent 
boson-boson loop diagram and the diagrams with three- 
body couplings are too weak to alter this behaviour. 

On the other hand, the main contributions for small cb 
come from the fermion and mixed fcrmion-boson loops, 
the latter arising in the three-body coupling. In partic- 
ular, the mixed boson-fermion loop diagrams containing 
the fermionic cut-off contribute to the evolution of the 
three-body coupling, even when the bosonic degrees of 
freedom have been integrated out. As a result, inclusion 
of A leads to a significant deviation from the mean-field 
result, as I <if = 2, that survives in the limit cb — > 0. 

These results for very large or very small values of cb 
should not be taken too seriously. Arguments based on 
"optimisation" of the cut-off function, see Ref. , indi- 
cate that one should choose the cut-off to try to maximise 
the rate of convergence for our expansion of the action. 
Although a precise criterion has not yet been defined for 
nonrelativistic theories, such arguments suggest a choice 
where bosons and fermion cut-offs run at roughly the 
same rate, i.e. cb is of the order of 1. 

In spite of this clear improvement over calculations 
that include two-body interactions only, adding the sim- 
plest possible three-body term is not enough to ensure 
that the results are completely independent of the pa- 
rameter cb in the region cb — 1. It is worth emphasis- 
ing again that, as long as any truncation of the effective 
action is made, the results will never be completely in- 
dependent of the choice of cutoff. It seems likely that 
further extensions of the effective action will result in 
stability of the results with respect to the variations of 
cs in wider region. Such extensions could include both 
four-body interactions as well as energy and/or momen- 
tum dependent three-body forces. Work along these lines 
is now in progress. 

In summary, we have performed an ERG analysis of 
the boson-boson scattering length in a system of non- 
relativistic fermions. Our study indicates that, while a 
simple ansatz with only local two-body interactions can 
yield results that are close to the exact value, these results 
are sensitive to the value of a parameter controlling the 
relative scales of the fermionic and bosonic cutoffs. Wc 
show that the inclusion of a local three-body interaction 
brings the scattering length closer to the exact value and 
significantly reduces its sensitivity to the relative scale 
parameter. 
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